Exact solution of close-packed dimers on the kagome lattice 
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It is well-known that exact enumerations of close-packed dimers can be carried out for two- 
dimensional lattices. While details of results are now known for most lattices, due to the unique 
nature of the lattice structure, there has been no complete analysis for the kagome lattice. Here 
we derive the close-form expression (1/3) \n(4:Xyz) for the free energy of close-packed dimers on the 
kagome lattice, where x, y, z are dimer weights. We use two different approaches: the Kasteleyn 
method of evaluating a Pfafhan and an alternative vertex model formulation. Both methods lead 
to the same final expression. The correlation function between two dimers at a distance equal or 
greater than two lattice spacings is found to vanish identically. 

PACS numbers: 05.50.+q,04.20.Jb,02.10.Ox 



I. INTRODUCTION 

A central problem in statistical physics is the enumer- 
ation of close-packed dimers on lattices. The origin of 
the problem has a long history tracing back to the 1937 
paper of Fowler and Rushbrooke [1] in their attempt of 
enumerating the absorption of diatomic molecules on a 
surface. A milestone in the history of the dimer prob- 
lem is the exact solution for the square lattice obtained 
by Kasteleyn [2] and Temperley and Fisher [3] in 1961. 
Indeed, the method of Kasteleyn is quite general and ap- 
plicable to all planar lattices [4]. Exact results obtained 
in this way are summarized in a recent review [5] for a 
number of two-dimensional lattices. 

In the case of the kagome lattice, however, there has 
been no complete analysis of the dimer problem other 
than studies of pure dimer enumerations most of which 
are numerical and series expansions (see [6] and refer- 
ences cited therein). In recent years there has been con- 
siderable interest in the study of physical phenomena on 
the kagome lattice. These range from high-T c supercon- 
ductivity [7], Heisenberg antiferromagnets [8-12], quan- 
tum dimers [13], to the occurrence of spin- liquid states 
[14]. It has also been shown that the consideration of 
close-packed dimers is related to the ground state of a 
quantum dimer model [15]. In light of these develop- 
ments, it is of pertinent interest to take a fresh look at 
close-packed dimers for the kagome lattice. 

In this paper we consider this problem, and derive the 
closed-form expression 
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(1) 



for the free energy (for definition of terms see below), a 
formula quoted in [5]. As exact solutions for other lat- 
tices are invariably of the form of a double integral akin 
to the Onsager solution of the Ising model [5], the very 
simple expression of the solution (1) and its logarithmic 
dependence on dimer weights are novel and unique. It 
points to the special role played by the kagome lattice 
(which often makes a problem more amenable), and sug- 



gests that caution must be taken in generalizing physical 
results derived from the kagome lattice. For example, 
as we shall see in Sec. IV below, the correlation between 
two dimers is identically zero at distances equal or greater 
than two lattice spacings on the kagome lattice, but this 
conclusion does not hold for other lattices. 

The kagome lattice is shown in Fig. 1, where x, y, and 
z are dimer weights along the three principal directions. 
We denote the lattice by C. Let A (= even) be the 
number of sites of C, so the lattice can be completely 
covered by A/2 dimers. The dimcr generating function 
is defined to be the summation 
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over all close-packed dimer configurations of C. Here, n x 
is the number of dimers with weight x, etc., subject to 
n x +n y +n z = A/2. Our goal is to evaluate the per-dimer 
free energy 
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in a close form. Past attempts have been confined to enu- 
merations of /(l, 1, 1). Here we consider the problem for 
general x, y, z, a consideration which can find application 
on anisotropic kagome system such as the volborthite an- 
tiferromagnet [12]. 

We derive the solution (1) using two different methods: 
the Kasteleyn method of evaluating a Pfaffian and alter- 
nately a method of a vertex model formulation, which wc 
describe in the next two sections. 



II. PFAFFIAN APPROACH 

The first step of the Pfaffian method is to find a Kaste- 
leyn orientation [2] of lattice edges. A Kasteleyn orien- 
tation of a planar lattice is an orientation of edges such 
that every transition cycle consisting of a loop of edges 
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FIG. 1: The kagome lattice C with dimer weights x, y, and z 
in the three principal directions. 




FIG. 2: The Kasteleyn orientation of the kagome lattice. A 
unit cell is the region bounded by dashed lines containing the 
six sites numbered 1, . . . , 6. 

derived from the superposition of two dimer coverings has 
an odd number of arrows pointing in the clockwise direc- 
tion, a property which we term as clockwise-odd. While 
Kasteleyn [4] has demonstrated that such an orientation 
is possible for all planar graphs, the actual orientation 
of edges for a given lattice, or graph, still needs to be 
worked out, and the crux of the matter of the Kasteleyn 
method is the finding of the appropriate clockwise-odd 
orientation. 

For the kagome lattice a Kasteleyn orientation can be 
taken as that shown in Fig. 2. The kagome lattice is com- 



posed of up-pointing and down-pointing triangles. The 
orientation in Fig. 2 consists of orienting all up-pointing 
triangles and every other down-pointing triangles in the 
counter-clockwise direction, with the other half of the 
down-pointing triangles oriented as shown. In this ori- 
entation a unit cell of the lattice consists of the 6 sites 
forming two neighboring down-pointing triangles, which 
are numbered 1, . . . , 6 as shown. Our orientation is dif- 
ferent from that used in [6]. 

To see that the orientation in Fig. 2 is indeed a Kaste- 
leyn orientation, we note that all transition cycles in 
Fig. 2 are clockwise-odd. As all transition cycles on the 
lattice can be formed by deforming those in Fig. 2 with- 
out altering the clockwise-odd property, all transition cy- 
cles are also clockwise-odd so the orientation in Fig. 2 is 
a good Kasteleyn orientation [16]. This is essentially the 
argument of Kasteleyn [4] . 

After rendering the Kasteleyn orientation, we can write 
the dimer partition function as a Pfaffian. As the analysis 
is now standard we follow the standard procedure [see, 
for example, Eq. (4.34) of [17]] and arrive at the result 
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lnPf[M(0,0)]d0d</> 
lndet|M(0,0)|d0d0 . (4) 



Here the factor 1/3 in the first line comes from the fact 
that there are 3 dimers per unit cell in a close-packed 
configuration, and Pf [M(6, 4>)] = ^dct|M(0, </>) is the 
Pfaffian of the 6x6 matrix 
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The a-matrices are read off from Fig. 2 to be 
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a(0,-l) = -a T (0,l) 

a(-l,-l) = -a T (l,l) 
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FIG. 3: An extended kagome lattice £ constructed by in- 
serting a decorating site attached to two inserted edges of 
weight 1 as shown. The decorating sites are denoted by solid 
circles. The unit cell is the region bounded by dashed lines. 
Repeating unit cells form a square lattice. 
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FIG. 4: Mapping between vertex and dimer configurations 
and the corresponding weights. 
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where the superscript T denotes the matrix transpose. 

The evaluation of the determinant in Eq. (4) gives the 
surprisingly simple result 



det[M(0,<£)] = 16x 2 y 2 z 2 



(7) 



The substitution of Eq. (7) into Eq. (4) now yields 
Eq. (1). The expression (7) and result (1) have previ- 
ously been obtained for x = y = z = 1 in [6] for pure 
dimer enumerations. 



III. VERTEX-MODEL APPROACH 

The kagome dimer problem can also be solved using 
a vertex model approach which is conceptually simpler. 
This involves the mapping of the dimer problem into a 
vertex model for which the solution is known. 

The first step of the mapping is to introduce the ex- 
tended kagome lattice £ of Fig. 3. The extended lattice 
£ is constructed from C by introducing 47V/3 extra lat- 
tice edges with weight 1 and 27V/3 new (decorating) sites 
as shown. By inspection it is clear that a bijection exists 
between dimer configurations on L and £ . This permits 
us to consider instead the dimer problem on £ . 

The dimer problem on £ is next mapped onto a vertex 
model. 



The extended lattice £ consists of N/3 unit cells each 
of which is the region bounded by dashed lines shown 
in Fig. 3. The unit cells form a square lattice S. We 
next map dimer configurations on £ into vertex config- 
urations on S, by regarding the 4 edges extending from 
a unit cell of £ as the 4 edges incident to a site on S. 
To each extending edge on £ covered by a dimer, draw a 
bond on the corresponding edge on S, and to each extend- 
ing edge not covered by a dimer, leave the corresponding 
edge empty. Then, as shown in Fig. 4, dimer coverings of 
a unit cell are mapped into vertex configurations on S. 
Since the number of bonds extending from each vertex is 
either 1 or 3, which is an odd number, we are led to the 
odd 8-vertex model considered in [18]. 

Vertex weights of the odd 8-vertex model can be read 
off from Fig. 4 as 



1*1 = XZ, 1*2 

u b = xy, i*6 
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(8) 



The per-vertex 8-vertex model free energy is then 
f 8v (x,y,z) = Jim — !— lnZ kagomo (x,?/, z) . (9) 

N^oc iV/o 

Comparing Eq. (9) with Eq. (3), we obtain the equiva- 
lence 



/kagomc(2!,y,20 = ( 3 ) fsv(x,y,z) . (10) 



Now the weights (8) satisfy the free-fermion condition 

UiU 2 + 1*3«4 = U 5 U 6 + U7l* 8 (11) 

for which the per-vertex 8-vertex model free energy has 
been evaluated in [18] as 

f8v= 16^J dd J d ^ lnF (^) ( 12 ) 
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The solution (1) is now obtained by substituting Eqs. (13) 
into Eqs. (12) and (10). 

IV. DIMER-DIMER CORRELATION 

The dimer-dimer correlation function can be evaluated 
by either considering a perturbation of the Pfaffian as in 
[19, 22] or by applying the Grassmannian method of [23, 
24], details of both approaches will be given elsewhere. 
Here we sketch steps in the Pfaffian computation. 

Define for the lattice edge connecting sites i and j in 
unit cell at r = (r x , r y ) an edge vacancy number 

nij.r = 1, if ij is empty, 

= 0, if ij is occupied, (14) 

where (•) denotes the configurational average. Then, the 
correlation function between two dimers on edges ij in 
cell ri and kt in cell r 2 is 

c(ij,r 1 ;k£,r 2 ) = {n, 3 . ri n k e, ri ) - («ij;ri)(nfc£;r 2 )- (15) 

To make use of Eqs. (15) we need to compute the dimcr 
generating function with specific edge(s) missing. Let A 
be the antisymmetric matrix derived from the Kasteleyn 
orientation, and let A' denote the antisymmetric matrix 
derived from A with edge ij, say in computing (n^), 
missing. Write 

Z = PiA 

Z' = PiA' = Pf [A + A] (16) 



where A is the matrix with zero elements everywhere 
except the ij element is —A.^ and the ji clement is 
-Aji (= Aij). Then 

(ny) = Z'/Z = Pi A' /Pi A (17) 

and 

. , 2 detA' det[A(I + G£)] , , T ^, AN 

™ = deTT - La = dct(/+GA) - (18) 

where G = A^ 1 is the Greens function matrix and / the 
identity matrix. 

In computing Eq. (18) we need only to keep those 
row(s) and column(s) in A and A' 1 where elements of 
A are nonzero. In addition, in the interior of a large 
lattice, the correlation depends only on the difference 
r = ri r 2 = {r x ,r y }, so elements of G are given by 

G(r) = J 2 J J^dOdcf, e^+^A-^e,^). 

(19) 

These considerations lead to the explicit evaluations of 
Eq. (18), and hence the correlation (15). 

Particularly, due to the fact that elements in A~ 1 (9, <f>) 
contain only a monomial of e ±l6 and e ±4< ^, a consequence 
of the fact that the determinant detA is given by the 
simple expression (7), the integral (19) vanishes for \ri x — 
T2x\ > 1 or \n y — r 2 y \ > 1. This leads to the result 

c(y,ri;fc*,r 2 )=0, |r 1 -r 2 |>2. (20) 

The absence of the dimer-dimer correlation beyond a 
certain distance, which is also found in the Suthcrland- 
Rokhsar-Kivelson state of a quantum dimer model [13], 
is a property unique to the kagome lattice. This under- 
scores the special role played by the kagome lattice in 
the statistical mechanics and quantum physics of lattice 
systems. 
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